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ABSTRACT 

Active Galactic Nuclei (AGN), Seyfert galaxies and quasars, are powered by 
luminous accretion and often accompanied by winds which are powerful enough 
to affect the AGN mass budget, and whose observational appearance bears an 
imprint of processes which are happening within the central parsec around the 
black hole (BH). One example of such a wind is the partially ionized gas respon- 
sible for X-ray and UV absorption ('warm absorbers'). Here we show that such 
gas will have a distinct signature when viewed in polarized X-rays. Observations 
of such polarization can test models for the geometry of the flow, and the gas 
responsible for launching and collimating it. We present calculations which show 
that the polarization depends on the hydrodynamics of the flow, the quantum 
mechanics of resonance line scattering and the transfer of polarized X-ray light 
in the highly ionized moving gas. The results emphasize the three dimensional 
nature of the wind for modeling spectra. We show that the polarization in the 
0.1-10 keV energy range is dominated by the effects of resonance lines. We pre- 
dict a 5 — 25% X-ray polarization signature of type-2 objects in this energy range. 
These results are general to flows which originate from a cold torus-like structure, 
located ~ lpc from the BH, which wraps the BH and is ultimately responsible 
for the apparent dichotomy between type 1 and type 2 AGNs. Such signals will 
be detectable by future dedicated X-ray polarimetry space missions, such as the 



NASA Gravity and Extreme Magnetism SMEX, GEMS (Swank et al. 2008). 



Subject headings: acceleration of particles - galaxies: active - hydrodynamics 
-methods: numerical - quasars: absorption lines - X-rays: galaxies 



Optical spectropolarimetry of Seyfert 2 galaxies was responsible for the revelation that 
both Seyfert 1 and Seyfert 2 AGNs are intrinsically the same objects, and that their apparent 
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differences are due to the inclination from which they are viewed ( Antonucci &: Miller|1985 ). 
A paradigm of a cold (T = 100 — 1000K and several parsecs in diameter) molecular torus 
which is obscuring the broad line region in most low and intermediate luminosity AGNs, is 
supported by the polarized optical line spectra and in some cases by direct interferometric 



observations (Jaffe et al. 2004). 



The obscuring torus may represent part of a flow originating at smaller (Elitzur & 



Shlosman 2006) or larger (Proga 2007) radii. Studies of the lines blueshifts combined with 



variability studies of warm absorber spectra (Behar et al. 2003; Netzer et al. 2003), suggest 



a location of the torus at the distance > lpc from the BH. 



First principal arguments (Krolik & Begelman 1986, 1988) support the idea that the 



torus is the source of an outflow which is evaporated by X-ray heating, and such gas is a 



plausible source of warm absorbers (Krolik & Kriss 1995, 2001). 



The question of the origin of the warm absorber flow is intrinsically related to the prob- 
lem of the vertical support of the torus against collapse. Infrared radiation pressure is a 



likely mechanism for such support (Krolik 2007), although magnetic fields could also con- 



tribute ( Lovelace et al.fl~998 ). One may also expect large scale magnetic fields playing a role 
in producing hydromagnetic flows related to warm absorbers. Models based on self-similar 



solutions of the type developed by Blandford & Payne (1982) were calculated by Emmer- 



ing et al (1991), Bottorff et al. (1997), and Bottorff et al. (2000); warm absorber spectra in 



transmission from a different type of 2D self-similar solutions of the Grad-Shafranov equation 



were calculated by Fukumura et al (2009). 



Warm absorber gas is heated and ionized by X-rays from the BH and in the low density 
limit such gas approaches the 'Compton temperature', T IC ~ 10 7 K, although multidimen- 
sional hydrodynamic simulations show that due to adiabatic losses the temperature may 
be smaller, ~ 10 6 K (Dorodnitsyn et al. 2008 a[b ). It also follows that mass loss rates 
of M = 0.1 — lM /yr and velocities of 100 — 1000 kms^ 1 are attainable, and that these 
properties match those of X-ray warm absorbers. 

Free electrons in the wind scatter radiation from the accretion disk and broad line 



region towards the observer giving rise to the polarized optical and UV (Antonucci & Miller 



1985) spectrum in Seyfert 2 objects. Highly ionized ions from this same mirror gas create 



rich X-ray spectra of warm absorbers at inclinations appropriate for Seyfert 1 objects. These 
correspond to photoionized plasma with numerous absorption and emission lines from highly 
ionized ions of Fe, Si, S, O, and Ne in the 0.1 - 10 keV energy range. These are detected 
by the Chandra and XMM — Newton satellites, in addition to features from many of the 
same ions in the UV seen by e.g. HST. Such spectra occur in a large fraction of low redshift 
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AGNs which are bright enough to allow detection (Reynolds 1997 McKernan et al. 2007). 



These lines occur in absorption in Seyfert 1 galaxies, and in reflection spectra in Seyfert 2 



galaxies (Kinkhabwala et al. 2002), corresponding to differing views of the flow and torus. 



X-ray polarimetry promises a new opportunity to study AGN winds through their X- 
ray reflection and thus to test AGN unification models. Here we report the results of self- 
consistent modeling of the X-ray polarization signature produced by a warm absorber flow. 
We choose a torus model with initial mass M tor = 9.3 • 1O 5 M , so the maximum Compton 
optical depth of the torus interior, r max = 40, is high enough to ensure capture of essentially 
all X-rays incident on the torus inner surface. We calculate the evolution of such gas and 
the fate of the torus, using the techniques described in Dorodnitsyn et al. ( 2008 a[b ), which 
includes the dynamics and the effects of photoionization and photon heating on the thermal 



and opacity properties in the gas (Kallman & Bautista 2001). 



Gas evaporated by X-ray heating is confined by the torus funnel. Our calculations suggest 
that to some extent, the throat of the torus acts similar to a de Laval nozzle in a rocket 
engine, simultaneously providing the flow with the fuel through evaporation of the material 
of the confining walls. This helps to convert the internal energy of the heated gas to bulk 
kinetic energy of the outflow. Some of the X-ray heated gas is close enough to the black hole 
to capture significant radiation flux and contributes to the polarization signature. This is 
enhanced by gas which fails to escape from the potential well, forming a returning current, 
as illustrated in Figure [TJ The opening of the torus funnel is affected by the hydrodynamics 
of evaporation and determines the fraction of AGN in which warm absorber spectra can be 



observed (Dorodnitsyn et al. 2008 b). 



2. Polarization signature of the X-ray excited flow 



Our calculations consist of three parts. The first part includes time-dependent, 2.5- 
dimensional hydrodynamic simulations, i.e. time-dependent 2-dimensional simulations with 



rotation. The equations of hydrodynamics are solved using the ZEUS-2D code (Stone & Nor 



man|1992 ), which we have modified to take into account various processes of radiation-matter 
interaction. The extensive heating occurs on a timescale which is much shorter than the dy- 
namical time, which makes the system of equations stiff, and thus including such processes 



requires a modification of the ZEUS hydrodynamical scheme (Dorodnitsyn et al. 2008 b). 



Radiation-matter interactions include: Compton and photoionization heating and Compton, 



radiative recombination, bremsstrahlung and line cooling. The XSTAR (Kallman & Bautista 



2001 ) code is used to calculate rates of heating-cooling, and also to calculate tables of line 
opacities. These are used in the final stage, when we calculate the emergent spectrum. We 
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Fig. 1. — The AGN torus, observed at t = 4 to (t is the dynamical time) from an inclination 
of 50°. Density is shown as the shaded equatorial ring and velocity streamlines are shown as 
colored tubes emerging from the central hole. Close to the axis the wind forms a fast bipolar 
outflow, with too low density and too high ionization to form observable warm absorber 
features. Favorable conditions for warm absorber viewing exist for inclinations 40 — 60°. 



-5 - 



calculate the transfer of polarized radiation in spectral lines in the rapidly moving medium 



adopting the Sobolev (1960) approximation, together with a three-dimensional version of 



the method of escape probabilities (Castor 1970 Rybicki & Hummer 1983). 



The gas evaporated from the surface of the torus by X-ray heating is confined within the 
torus funnel within a vertical range of z ~ 1 — 2 pc. The sonic surface is located at, or very 



close to the surface of the funnel (see Fig. 3 of Dorodnitsyn et al. (2008 b), where the flow 
dynamics is described in detail). The torus maintains aspect ratio of R/H ~ 1 until very late 
times when it is significantly depleted by the mass loss. Simulations show that the flow is very 
highly ionized close to the axis with a temperature approaching Compton temperature, and 
at higher inclinations as the wind is more dense, self-screening provides favorable conditions 
for the formation of warm absorber spectra. At inclinations where warm absorber gas is 
observed, the line-of-sight first touches the stripped cool gas of the torus outskirts, and then 
penetrates through the torus wall. The distribution of the ionization parameter, £ varies 
dramatically along the line-of-sight. At inclinations, 9 > 50° — 60°, the radiation of the 
nucleus is heavily obscured. Warm absorber spectra are produced within a range of angles 
A6 ~ 10° — 15° at such inclinations. The inferred velocity of the gas decreases as the 
inclination increases. 

The results of hydrodynamic modeling are post-processed with our radiation transfer 
code. This code makes use of the Sobolev approximation to calculate a three-dimensional 
transfer of the polarized radiation in approximately 10 5 spectral lines in 0.1 -10 keV range. 
A precursor of this approach included 3D Sobolev radiation transfer of the un-polarized 



radiation and is described in Dorodnitsyn & Kallman (2009). Scattering of the polarized 



radiation in a resonant spectral line with corresponding calculations of the scattering phase 



matrix, was studied by, e.g. Chandrasekhar 



( 


1960 


)• 


Hamilton 


( 


1947 


), 


Lee et al. 


( 


1994 



However, calculations which take into account transfer of polarized light through a general 
3-dimensional distribution of opacity are numerically demanding. To make them tractable 
for the many lines in a real gas, significant modification is required both to the methods and 
the algorithm. 

Following Chandrasekhar, we describe the intensity of the polarized radiation by the 
Stokes vector: I = (Jl, Jr, U, V), where Jl, Ir are the intensity components projected on 
two arbitrary, orthogonal axes; U measures the polarization 45° away from the axis, and 
V measures circular polarization. The equation, describing the transfer of the polarized 



radiation in arbitrary direction s, reads, (e.g. Chandrasekhar (I960)): 



n 



°°,r,n) 



tfl(i/°°,r,n) 
ds 



[XW) (S(z/ 00 ,r,0)-I(z/ 00 ,r,0)), 
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where S is a Stokes vector source function, v°° is a frequency of the photon, measured in the 
"lab" frame, and (p is the spectral line profile, and x\ is the line opacity. Every resonance 
with a spectral line, is accompanied by a corresponding change in the radiation field: 



Iafter = Ibefor e e TM + S(l - C ^) = I d ist,i + I]oc,i, (2) 

that is, the radiation field consists of the local and distant contributions; the Sobolev optical 



depth r\ (Castor 1970; Rybicki & Hummer 1978 Dorodnitsyn & Kallman 2009) is assumed 



to be independent of polarization. The first term in (|2j gives the attenuation of the external 
radiation and the S(l — e~ ru ) term is the reinforcement of the radiation due to scattering 
of the radiation from all directions into the direction towards an observer. 

The Sobolev optical depth scales the escape probability, /3 n which is given by 

P n = (1 _ e -M(n))/ e -u(n)_ (3) 

Integrating we adopt the disconnected approximation (Grachev & Grinin 1975; Marti & 



Noerdlinger||1977| Rybicki & Hummer 1978, 1983); that is, we assume that the interaction 



between resonances along the line-of-sight is happening only in the forward direction (towards 
the observer), and that all other possible interactions are neglected. Only resonant line 
transitions are considered, and complete frequency redistribution in the co-moving frame is 
assumed. 



The source function, S is calculated from the following (Hummer 1969): 



S(isn) 



47T 



R(z/n', i/n) J(i/n) di/dn', 



(4) 



where R(i/n', un) is the scattering phase matrix, given by |Chandrasekhar~ (1960). This 
matrix is a generalization of the case of a linear combination of the Rayleigh and isotropic 
phase matrices; coefficients of the matrix depend on the angular momentum quantum num- 
bers of the initial and final states of the transition (Hamilton 1947). The mean intensity, 



J(z/) depends on the availability of radiation, and thus on a source function, which itself 
depends on angles through its dependence upon /3 n . 

We take approximately into account the depolarizing effect of finite Sobolev optical 
depth, and thus avoid overprediction of the polarization fraction. That is, when the local 
escape probability, /3 esc << 1, the photon suffers many uncorrelated resonance scatterings 
within the resonant region, which tend to depolarize and make isotropic the local contribution 
to the radiation field. Solving such a problem in a medium with spatial velocity gradients is 
better done with a local Monte Carlo simulation. However, for our purposes, that would be 



-7- 



numerically prohibitive. To make the problem tractable, we adopt a simplifying approach in 
which the source function is represented as a combination of the source function derived for 
an isotropic phase matrix for the local contribution to the radiation field, and an anisotropic 
part for the distant contribution. The mean intensity J consists of two parts: one is the local 
contribution which depends on the source function, S, and the other is given by the radiation 
coming from the external source (the "core"). The idea is to calculate J, form the angular 
independent, un-polarized source function, , and then to substitute it in (pi). Such an 
approximation captures important properties of Sobolev transfer of the polarized radiation, 
anticipated also on physical grounds. Making use of the method of escape probabilities 



(Castor 1970), S is found with the further help of the relation: J = S^°^(l— < /3 n >4 n )+ < 
I c /3 n >n SIC , where we make use of the notation: < / >^= (47r)~ 1 j Q /(f2') dQ! . The last term 
represents the probability for the radiation of the core, I c to reach a given point, weighted 
by the probability of penetration. We assume that the core emits continuum radiation of 
specific intensity J c which is constant near the frequency of the line. For the source function 
which takes into account the effects of polarization we obtain: 

S«(n) = J R(n» {S<°>(1 - M+ < IcA' >n SIC } dtt, (5) 
where the un-polarized source function, So is found from the relation: 

c(o) _ 1 < IgPn >n BIC / R x 

In case of isotropic scattering, from ^ one will recover (|6]), i.e. the un-polarized source 
function, which is usually used in Sobolev calculations of the line spectra. The other extreme 
is (3 n > = 1; all photons arriving from the core are scattered and leave the resonance region 
without any further interaction, so that no gradient-related effects are important. On the 
other hand, if (3 n i << 1 a given photon is scattered many times within the same local 
resonant region. The resultant intensity may get significantly depolarized. 

The resultant polarization of the local contribution to the intensity is largely determined 
by the last scattering event, i.e. by the R(n', n) term. Further details of the method 



including the geometry are given in Dorodnitsyn & Kallman (2009) for the non-polarized 



case (for the geometry of the problem compare with their Figure 1). Here we formally 
generalize to the case where four components of the Stokes vector are transfered. We have 
rigorously tested our code with respect to problems of Sobolev transfer of the polarized flux 
in individual lines and the distribution of lines. Prescribing different distributions of the 
velocity and density we are able to compare transfer of the individual polarized modes with 
results of un-polarized transfer calculations. 
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Fig. 2. — Warm absorber spectra viewed at inclination 50° in two orthogonal linear polar- 
ization modes, IlJr (up). Polarization fraction, P (bottom). 
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Fig. 3. — Same model as in Figure [2] but at inclination 70°. 
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3. Results 

A notable feature of most warm absorber spectra is the fact few of the lines appear 
truly black in their cores. While this could be due to finite resolution or finite column of 
the gas, in our simulations this is because of scattering into the line core. Such scattering 
typically occurs at angles close to 90°, and so the residual flux in the lines is highly polar- 
ized. This is similar to the polarization of the (rest frame) UV lines observed from broad 



absorption line quasars (Ogle et al. 1999). Figure |2| shows spectra and the corresponding 
polarization fraction, P at inclination 9 = 50° after the wind simulation has evolved through 
four dynamical times (the dynamical time is given by: to — 1-5 x 10 4 r pc (M-qu/Mq)~ 1 I 2 (yr), 
where M-qk/Mq is the BH mass in units of 10 6 M Q and r pc is the radius of the torus in pc) 
Each synthetic spectrum is convolved with a Gaussian with width and centroid related by 
<je — 10 -3 .E. In the 0.1-10 keV range, the most prominent feature of the spectrum is an 
absorption trough (Figure |2j upper panel) and a corresponding maximum in polarization at 
~ 14.9 — 17.2A, (0.72 — 0.83 keV). This is due to H-, and He-like oxygen and moderately ion- 
ized iron, FeXVII — FeXIV. The polarization maximum has P max = 10%, and polarization 
at FWHM=0.05 keV about 2%. This is produced at a distance ~2.1 pc from the BH, where 
the radial column density rises from iV co i = 1.5 • 10 22 to iV col = 3.8 • 10 22 , and the ionization 
parameter, £ ( Tarter et al.||l969 ) drops from 6.4 • 10 3 to 500 — 10 within a region 2.8 — 4.4 pc. 



This results in an ideal situation for scattering in the flow, since further away the plasma is 
again over-ionized: £ ~ 10 4 . Polarization of lines in the 0.1-2 keV band is generally <1% 



0, 



but for the strongest lines the polarization can reach 2 - 3% in the O VIII Let and Ne IX 
He a lines; and ~10% in the case of the features at ~ 0.75 keV. 

At high inclinations, 9 > 60°, most lines appear in emission, and the continuum flux 
is ~ 10 -2 of that at low inclinations. Figure [3] shows the spectrum and corresponding 
polarization fraction, P at inclination 9 = 70°. One of the insights provided by multi- 
dimensional transfer calculations is the importance of blending of multiple lines, i.e. when 
several lines collectively appear as a single strong feature. Blending of multiple lines from Fe 
XXV and Fe XXVI forms such a spectral feature at Eq ~ 6.5 keV, at the inclination in figure 
[3j resembling a P-Cygni line. However, this shape is not due to the effect of the moving wind, 
rather it is an interplay between geometric effects and ionization balance. That is, due to the 
attenuation of the continuum at high inclinations, when the continuum intensity in the line 
emitting region falls below a threshold, I c < S\(l — e T ), where Si is the source function and 
T\ is the optical depth, then the ensemble of lines corresponding to lower ionized gas pops 
up above the continuum, while the more highly ionized gas is still seen in absorption. Other 
examples of such features occur at lower inclinations, such as at 9 ~ 45°, corresponding to 
the O VII edge at 16. 8A (739 eV), together with lines from the L shell of Fe. 
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Figure [4] shows the distribution of polarization fraction versus inclination angle. We 
plot the maximum polarization, P max (P convolved with a Gaussian with oe = (15)~ l E), 
and the flux mean < P >p in the 0.1-2 and 2-10 keV bands. This shows that the maximum 
polarization is < P > max ~ 15 — 20% over a wide range of angles. The various peaks and 
troughs in this plot are reflections of the density structure in the flow, including wedges 
of gas which are stripped off the torus. The flux mean polarization has a maximum near 
9 = 67°, where < P > F (0.1 - 2) keV < 10% and < P > F (2 - 10) keV < 2%. At 9 > 75° 
the flux mean polarization in the 0.1-10 keV range is < P >f^$ 20% and the maximum is 
-Pmax — 70%. At such inclinations the torus is Compton thick along the line of sight and the 
continuum flux is very low, < 1% of the unattenuated value. 

Scattering by free electrons from optically thin plasma can produce a continuum polar- 

. Optically thin, 

single scattering approximation provides an upper limit for the polarization. This approach 
takes into account the geometry of the problem together with optical depth effects and pre- 
dicts both the polarized flux and the fractional polarization is less than those produced by 
lines. This is shown as the labeled curve in Figure [4]). At those viewing angles where the 
Thomson continuum has significant polarization, the flux is ~ (10~ 3 — 10 _4 )F C , where F c is 
the flux of the core. 

Our results show that models in which the obscuring torus is the source of an X-ray 
excited flow generate potentially observable polarization signatures due to significant scat- 
tering of X-rays by atomic features in the flow. Thomson electron scattering continuum at 
9 < 60° has polarization < 1%; it is increasing to a maximum of 17% at 9 ~ 80°, and then 
decreases at larger 9. Polarization from lines is greater and has the following properties: at 
very low inclination, 9 < 45°, corresponding to prototypical Seyfert 1 galaxies, the linear po- 
larization fraction is negligible averaged over the spectrum, but most absorption lines have 
some residual, highly polarized flux; the maximum polarization fraction in strong lines is 
1 — 10%. At higher inclinations, 45° < 9 < 75° in the 0.1-2 keV range the mean polarization 
is 10%, and the maximum polarization is 20-30%; in the 6-8 keV range the mean polarization 
is 2-10%, and the maximum polarization 20-30%. At still higher inclinations 9 > 70 — 80° our 
models are very Compton thick. Scattering in X-ray lines above the torus funnel produces 
very strongly polarized flux with mean polarization fraction more than 20% and maximum 
polarization 60-70%. Although these predictions are specific to our dynamical calculations, 
the general behavior will hold for any scenario involving an outflow which is shadowed by a 
torus with aspect ratio near unity. Observations which fail to show such high polarization 
would call into question current ideas about the torus location, shape and thickness. 



ization fraction which is comparable to that in the lines (Bianchi et al.||2009 
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Fig. 4. — Distribution of the polarization fraction with inclination angle. Curves: dashed: 
maximum polarization fraction within the energy range; solid: thin: flux mean polarization 
fraction within the energy range; thick: Thomson scattering continuum 
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